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INTRODUCTION 

In  this  report  we  present  a  shock  spectrum  approach  to  the  design  of  lightly  damped  rela¬ 
tively  light  equipment  in  structures  subjected  to  ground  shock.  The  analysis  uses  as  a  model  an 
jV-degree-of-freedom  structure  to  which  is  attached  a  single-degree-of-freedom  component.  It 
is  known  that  significant  interaction  effects  occur  when  the  equipment  frequency  is  close  or 
equal  to  one  of  the  natural  frequencies  of  the  structure,  referred  to  as  tuning.  Otherwise,  the 
conventional  floor  spectrum  method  which  neglects  interaction  is  valid. 

In  previous  research  the  undamped  equipment-structure  problem  for  tuned  systems  [1] 
and  slightly  nontuned  systems  [2]  was  studied  and  it  was  shown  that  a  design  ground  response 
spectrum  together  with  fixed-base  dynamic  properties  of  the  structure  alone  and  of  the  equip¬ 
ment  alone  can  be  used  to  estimate  the  peak  response  of  the  equipment.  The  analysis  took 
advantage  of  the  mathematical  structure  of  the  equations  and  of  asymptotic  methods  made  pos¬ 
sible  by  the  smallness  of  the  equipment  mass  in  comparison  with  the  mass  of  the  structure  to 
obtain  simple  results  for  tuned  and  nearly  tuned  systems.  In  this  paper  we  extend  the  results  to 
damped  tuned,  nearly  tuned,  and  completely  nontuned  systems.  For  light  equipment  and  small 
damping,  the  results  obtained  can  be  easily  and  efficiently  implemented  by  a  designer.  The 
most  important  aspect  of  the  analysis  is  its  extreme  simplicity;  namely,  if  the  response  spectrum 
for  the  ground  motion  is  available,  the  peak  response  of  the  equipment  can  be  calculated 
merely  by  multiplying  the  former  by  an  amplification  factor. 

This  approach  is  in  contrast  to  several  earlier  analyses  of  equipment  response,  such  as  the 
floor  spectrum  method  in  which  the  equipment  is  treated  as  a  single-degree-of-freedom  system 
subjected  to  a  base  motion  that  is  taken  to  be  that  which  the  structure  would  experience  at  the 
attachment  point  in  the  absence  of  equipment.  Not  only  does  this  method  neglect  interaction, 
it  has  the  further  disadvantage  of  requiring  that  an  expensive  time-history  analysis  of  the  struc¬ 
ture  be  conducted  in  order  to  determine  a  base  motion  [3].  Approximate  techniques  that 
bypass  associated  computational  problems  have  been  proposed  whereby  floor  response  spectra 
are  developed  from  ground  spectra,  but  these  methods  are  ad  hoc  and  their  accuracy  cannot  be 
evaluated. 

It  is  also  possible  to  account  for  interaction  by  considering  the  system  to  be  an  N+l- 
degree-of-freedom  which  is  then  treated  by  a  modal  analysis  using  an  ensemble  of  spectrum- 
consistent  ground  motion  records.  Neither  approach  is,  however,  ideal.  If  equipment-structure 
interaction  is  important,  two  closely  spaced  modes  will  appear,  the  contributions  of  which  must 
be  summed  when  a  modal  approach  is  used,  but  there  exists  no  consensus  as  to  an  appropriate 
summation  procedure.  A  disadvantage  of  the  second  approach-time-history  analysis-is  that  it 
is  very  expensive  and  for  a  given  item  of  equipment  must  be  carried  out  for  a  wide  range  of 
ground  motion.  If  the  equipment  is  light,  the  use  of  an  iV+l-degree-of-freedom  model  with 
standard  structural  dynamics  computer  codes  may  mask  significant  response.  The  method 
developed  here  allows  the  design  engineer  to  use  the  given  ground  motion  design  spectrum 
directly  to  determine  the  response  of  the  equipment.  The  expense  of  time-history  analysis  and 
the  difficulties  of  evaluating  response  when  closely  spaced  modes  are  involved  are  thus  avoided. 

In  this  report  we  describe  physical  and  numerical  experiments  on  the  response  of  tuned 
and  nearly  tuned  equipment-structure  systems.  The  physical  experiments  were  carried  out  on 
the  large  shaking  table  at  the  University  of  California,  Berkeley,  and  involved  a  40,000  lb, 
three-story  steel  frame  model  with  oscillators  simulating  items  of  equipment.  The  model  was 
subjected  to  a  wide  variety  of  simulated  ground  motions  and  the  peak  accelerations  in  the 
equipment  items  measured.  The  experiments  confirmed  the  accuracy  of  the  analysis  in  predict¬ 
ing  the  response  of  light  equipment. 

In  addition  to  the  physical  experiments  we  have  performed  numerical  experiments  using 
as  a  model  a  standard  ten-story  building  and  a  single  equipment  item.  The  numerical  experi¬ 
ments  corroborated  the  results  obtained  in  the  physical  experiments  and  indicated  some  surpris¬ 
ing  aspects  of  conservatism  in  generally  accepted  modal  summation  rules.  It  is  shown  in  this 
report  that  both  the  SRSS  summation  method  and  the  NRL  summation  method  can 
significantly  overestimate  equipment  response  in  damped  tuned  systems. 
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The  results  of  the  analysis  for  the  lightly  damped  light  equipment  in  ground  shock  loaded 
structures  presented  here  can  be  used  to  design  the  damping  in  the  structure  and  in  the  equip¬ 
ment  to  minimize  the  response  of  the  equipment.  It  is  shown  that  the  best  choice  of  damping 
is  when  the  damping  in  the  structure  and  the  equipment  are  equal. 

The  results  developed  here  have  been  mainly  concerned  with  short  duration  loading  typi¬ 
cal  of  ground  shock,  but  by  application  of  random  vibration  theory  to  the  tuned  equipment- 
structure  system  it  can  be  shown  that  for  very  long  duration  ground  motion,  similar  predictions 
of  equipment  response  can  be  developed.  These  results  could  be  useful  in  predicting  the 
response  of  light  equipment  in  large  ground  vehicles  traversing  randomly  rough  track  such  as 
the  MX  system. 


ANALYSIS 

The  equations  of  motion  of  the  N-degree-of-freedom  structural  system  to  which  is 
attached  a  single-degree-of-freedom  equipment  item  (Figure  1)  take  the  form 

(Mij  Uj  +  Cij Uj  +  K,j Uj)  -  jt  (CjjRjiig  +  KjjRjUg)  +  Fe,  ,  i-\,2,...,N  ( 1 ) 

j->  J- i 

where  M0,  C,j,  and  Ku  are  the  mass,  damping,  and  stiffness  matrices,  respectively,  and  U,  is 
the  absolute  displacement  of  the  i'h  degree  of  freedom.  The  vector  R,  is  a  vector  of  influence 
coefficients  introduced  to  couple  the  ground  motion  ug(t)  to  the  structure,  and  e,  a  vector 
whose  components  are  zero  at  every  degree  of  freedom  except  the  one  to  which  the  equipment 
is  attached,  denoted  by  the  index  r,  where  it  takes  unit  value.  The  term  F  is  the  interaction 
force  between  the  equipment  item  and  the  structure.  The  equation  of  motion  for  the  equip¬ 
ment  displacement  u  is 

-mu  -  F  =  c(u-Ur)  +  k(u~Ur)  (2) 


where  m,  c,  and  k  are  the  mass,  damping,  and  stiffness  of  the  equipment,  respectively.  Using 
classical  modal  description  of  the  structure  and  Laplace  transform  techniques,  the  solution  for  u 
for  the  multidegree-of-freedom  system  takes  the  form 


( p 2  +  2 p<ap  +  <u2)  +  p2  J) 
*-i 


/nd>*2(2)3to/7+to2) 
Mk(p2+2Bk(lkp+ill ) 


4>  *  Z  ♦,*  Z  M„R,(2Bk  n  kp+Sl  k2)  (2/3 u>p+<o2) 

/-I  /-i _ 

Mk(p2+2BkSlkp+(lf) 


(3) 


where  /3  is  the  fraction  of  critical  damping  for  the  equipment,  n*  is  the  k'h  natural  frequency  of 
the  structure  alone,  4>*  is  the  i'h  component  of  the  k,h  mode  shape  of  the  building  alone,  Mk  is 
the  generalized  mass  for  the  k'h  mode,  and  Bk  is  the  fraction  of  critical  damping  in  the  k'h 
mode  assumed  small  enough  not  to  couple  the  modal  equations. 

The  solution  of  eq.  (1)  will  be  obtained  by  residue  theory.  The  poles  of  the  transfer  by 
which  u  is  related  to  ug  are  simple  poles  and  determined  by  the  zeroes  of  the  denominator  of 
the  expression  in  eq.  (1).  Two  different  situations  arise,  one  when  the  equipment  is  tunecl  or 
nearly  tuned  and  the  other  when  the  natural  frequency  of  the  equipment  is  well  away  from  any 
of  the  structural  frequencies  i.e.,  is  grossly  nontuned.  In  both  cases  the  poles,  for  light  equip¬ 
ment,  appear  near  the  natural  frequencies  of  the  structure  alone  and  of  the  equipment  alone. 
For  tuned  or  nearly  tuned  equipment,  two  closely  spaced  poles  (tuning  poles)  appear,  located 
near  the  equipment  frequency  and  the  frequency  of  the  structure  to  which  the  equipment  is 
nearly  tuned,  one  below  those  frequencies  and  one  above  them.  For  an  undamped  perfectly 
tuned  system,  these  two  poles  coalesce  into  a  double  pole  as  the  equipment  mass  approaches 
zero.  This  gives  rise  to  a  response  in  the  time  domain  that  grows  without  bound  proportional 
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to  time  regardless  of  the  time  history  of  the  input.  Thus,  the  contribution  to  the  sum  of  the 
residues  at  all  poles  is  dominated  by  the  residues  at  the  tuning  poles.  The  contribution  of  the 
summation  term  to  the  residues  at  these  two  poles  is  dominated  by  the  term  k  =  n  since  the 
denominator  of  that  term  is  nearly  zero.  Hence,  in  the  region  of  the  tuning  poles,  which  is 
around  p  =  io>,  eq.  (3)  can  be  approximated  by 


U 


(p2  +  2fia>p  4-  a>2)  +p2y 


2f}<i>p+(a 2 

p2+2BnCl„p+(l2 


( 2/3ajp+a>2)(2B„nnp+n  2)  ^  „ 

p2+2B„nnP+a2 


(4) 


where 


y _ * -  (5) 

MJ  2 

and 

c;  =  *;£  £  (6) 

are  the  effective  mass  ratio  and  participation  factors,  respectively. 

The  contribution  to  the  response  of  the  nontuned  poles  is  quite  standard.  Their  locations 
are  close  to  what  they  would  be  for  the  structure  alone,  namely 

p  =  ~Bmtlm  ±  inm  (7) 


Evaluating  the  residues  at  these  poles  and  dropping  negligible  terms,  we  obtain  to  dominant 
order  the  contributions  from  the  nontuned  structure  poles: 

t  /v 

£,  n  n  /  nme~fimn,"(^r)sinft„,(f-T)t/T  (8) 

0  m-1 

m^n 


The  contribution  to  the  response  from  the  tuning  poles  is  different  and  must  be  obtained  by  a 
special  approach.  This  has  two  parts;  a  dominant  portion  that  will  be  developed  subsequently 
and  that  has  a  later-occurring  peak,  and  a  nondominant  portion  that  is  of  the  same  order  as  the 
terms  above.  The  second  portion  is  easily  obtained  and  takes  the  form 

'r  N  Cm 

I  uAt)  £  — - 7 - —  caf>'eu'('_T)sin w(t~T)dr  (9) 

*b  *  ",  (l-w/flj2 

nr&n 


The  dominant  portion  of  the  contribution  from  the  tuning  poles  given  by  eq.  (4)  may  be  writ¬ 
ten  in  the  form 

u  =  [N(p)/D{p)]ug  (10) 

where 

Nip)  -  (2/3wp+w2)l2fi„(l+^)cup  +  (l+f):w2]C,n  (11) 


Dip)  -p4  +  <up,l2/3(l+y)+2B„(l+?» 

+  wyi2  +  y  +  2*  +  £2-M0fl,,(l-K)] 

+  to3p  (2/3  ( 1  ) 2  -h  2  (!+£)]  +  to4(l+^)2  (12) 
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In  the  above,  £  -  (Q„-ai)/co  is  the  nontuning  parameter. 

The  roots  of  Dip)  will  be  close  to  ±  iw  since  y,  /3,  B„,  and  (  are  taken  to  be  small. 
To  locate  the  poles  of  D(p),  we  replace  pin  Eq.  (12)  by 

/>-/»(  1+8)  (13) 

where  8  is  a  small  quantity  in  terms  of  which  eq.  (12)  becomes 
84  +  {4—  /[2/3(l+y)  +  25„(l+£)]}83 
+  {4-y-2{-£2-40fln(l+£)-il6/3(l+y)+6B„(l+£)lU2 
+  {— 2y - 4(  - 2f2 -  8/8BJ1+0  -  / [2/3(2+3y-2£-£J)  +  4£„  (1+f )]}  8 
+  l-y-4/33„(l+e)  ~  'l2/3(y~2f-e2)1)  =  0  (14) 


The  solution  of  Eq.  (14),  to  dominant  terms,  is 

s  =  f  ±y+/(^T^±f)  (15) 

where  here  and  throughout  the  remainder  of  the  analysis  the  upper  signs  are  taken  together  to 
give  one  root  and  the  lower  the  other.  The  quantities  X  and  p  are  given  by 

X  =  ^[{[y  +  e-(fi-ByY  +  4e(lS-Bn)2Y/>+[y  +  e-(ii-Bn)2]\:  (16) 

P  -  ^({(r  +  |2-(i3-5fl)J]2  +  4^(i3-fin)2),/'-[y  +  e2-()8-Bn)2]p  (17) 


The  formal  inversion  of  (10)  is 


Hit)  - 


_1_  f 

2 ni  { 


mei 

Dip) 


uAp)ep'dp 


(18) 


where  T  is  a  suitable  Bromwich  path.  If  ugip)  is  taken  to  be  1,  then  the  inversion  directly 
yields  that  portion  of  Green’s  function  ucit)  for  the  solution  which  comes  from  the  dominant 
contribution  of  the  tuning  poles.  The  solution  for  the  equipment  acceleration  for  given  ground 
motion  Hgit)  from  the  dominant  contributions  of  the  tuning  poles  takes  the  form 


uit)  «  f  Hcit-^H^Mdr 


(19) 


The  function  HG  will  be  obtained  by  residue  theory,  since  there  are  no  branch  cuts  in  the  p 
plane.  The  inversion  of  uc  for  the  general  case  (eq.  (15))  is  obtained  by  writing  the  denomina¬ 
tor  Dip)  in  the  form 

Dip)  -  ip~p])ip-p])ip-p7)ip-p2) 


where 

p i  “  (1  +  j  +  y2 —  w  ' — 2 —  +  2  ' 

.  X,  ,P+B„  u. 

p2  ”  /w  (1  +  j  —  y)  ~  w  j - 2 

and  p\  and  p2  are  the  complex  conjugates  of  p{  and  p2.  Evaluating  the  residues  at  each  pole 
and  collecting  complex  conjugate  terms  in  pairs,  we  obtain  the  result,  correct  to  dominant 
order, 
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wc0)  =  C"  ,  e  <^+fil"'/2  [X  sinh-^w/ cos-^-wf  sin(l  + 

2  2  2 

-  X  cosh-^<o/  sinytuf  cos(l  4-  y  )a>f 

-  n  sinh^to/ cosy d)t COS (1  +  y)w/ 

-  p.  cosh^-oi/  sinytu/sin(l  +  y)o>f)  (20) 

The  dominant  contribution  of  the  tuning  poles  to  the  acceleration  response  u(t)  to  a  specified 
imposed  ground  motion  ug(t)  is  obtained  by  substituting  the  above  equation  in  Eq.  (19).  The 
complete  solution  for  the  response  of  the  equipment  is  then  obtained  by  summing  eqs.  (8), 
(9),  and  (19)  and  takes  the  form 


uCO-fu/r)  £  r- 
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cue  T,sin&>0— t)  +  wc(r— t)I</t 


where  uG  is  given  by  eq.  (20). 

The  character  of  the  two  parts  of  the  solution  in  Eq.  (21)  differs.  The  contributions  from 
the  nontuning  poles  and  the  nondominant  contributions  from  the  tuning  poles  are  conventional 
and  would  attain  their  peaks  during  the  ground  excitation  or  shortly  thereafter.  The  dominant 
response  from  the  tuning  poles,  on  the  other  hand,  is  controlled  by  the  energy  transfer  from 
the  structure  to  the  equipment  through  beating,  which  takes  a  relatively  long  time  [1,  2],  and 
the  peak  response  of  this  contribution  will  occur  substantially  after  the  end  of  the  ground 
motion. 

The  peak  response  of  grossly  nontuned  systems,  i.e.,  where  the  equipment  frequency  is 
far  from  all  structural  frequencies,  can  be  estimated  in  the  same  way.  For  light  equipment,  the 
structure  poles  are  only  slightly  shifted  from  their  location  for  the  structure  alone  and  additional 
poles  due  to  the  equipment  occur  close  to  those  for  the  equipment  alone: 

p  =  — /3ai  ±  iw  (22) 

The  residues  at  the  structure  poles  are  as  before.  The  residues  at  the  equipment  poles  contri¬ 
bute  a  term  similar  to  the  nondominant  term  in  the  tuned  case.  The  derivation  is  standard  and 
the  complete  response  for  the  equipment  in  the  grossly  nontuned  case  is  given  by 
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APPROXIMATE  RESPONSE  SPECTRA  RESULTS 

The  results  given  in  the  previous  section  can  be  applied  to  the  design  of  equipment  and 
equipment  mounting.  Typically,  the  information  available  to  the  designer  is  a  ground  shock 
spectrum.  It  was  shown  in  references  1  and  2  that  for  the  case  of  the  undamped  system,  both 
tuned  and  nearly  tuned,  the  peak  acceleration  in  the  equipment  could  be  related  to  the  response 
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ft 


spectrum  of  the  ground  motion.  The  present  paper  is  concerned  with  the  generalization  of 
these  results  to  the  damped  tuned  or  nontuned  cases, 


For  the  grossly  nontuned  system  where  the  equipment  frequency  is  far  away  from  all 
structural  frequencies  and  these  are  well  separated,  conventional  summation  methods,  e.g.,  the 
square  root  of  the  sum  of  the  squares  procedure  applied  to  eq.  (23),  lead  to  an  estimate  for 
Hmax  'n  the  form 
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(24) 


where  SA(fl, B)  is  the  acceleration  ground  shock  spectrum  for  the  ground  motion  uK  evaluated 
at  frequency  fl  and  damping  factor  B.  Other  responses  can  be  estimated  by  substituting  the 
appropriate  ground  shock  spedtra  in  the  above  result. 

The  approach  for  tuned-undamped  systems  is  described  in  detail  in  reference  1.  This  lejd 
to  the  following  estimates  of  maximum  acceleration  and  displacement 


u 


Imax 


S^(w,0) 


(25) 


and 


Sd{io,  0) 

y'h 


(26) 


where  SA((i>,p),  SD(a>,p)  are  the  acceleration  and  displacement  ground  shock  spectra  for  fre¬ 
quency  «j  and  damping  factor  p. 

The  simplicity  of  the  result  can  be  explained  on  physical  grounds.  In  weakly  coupled  sys¬ 
tems  with  the  same  frequency,  the  response  of  the  system  involves  a  perfect  energy  exchange 
between  each  component  at  a  beat  frequency  much  lower  than  the  natural  frequency  of  each 
component.  The  same  phenomenon-a  classical  beat  phenomenon-occurs  here.  The  coupling 
is  weak  because  the  ratio  of  equipment  mass  to  structure  mass  is  small. 

When  a  structure  is  subjected  to  a  ground  motion,  the  velocity  imparted  to  the  structure 
is  mass  independent  and  determined  only  by  the  ground  motion.  Thus,  if  the  same  ground 
motion  were  applied  directly  to  tuned  equipment,  the  same  velocity  would  be  transmitted  to  it. 
Kinetic  energy,  on  the  other  hand,  is  proportional  to  the  mass  of  the  system  excited;  in  equip¬ 
ment,  that  energy  would  be  much  smaller  than  in  a  structure.  However,  if  the  equipment  were 
attached  to  a  structure  and  the  structure  subjected  to  a  ground  motion,  the  kinetic  energy 
imparted  to  the  latter  would  be  wholly  transmitted  to  the  equipment,  if  tuned,  and  the  velocity 
imparted  would  be  amplified  by  the  reciprocal  of  the  square  root  of  the  mass  ratio.  Damping  is 
clearly  important  in  this  process  because  the  energy  transfer  requires  many  cycles  and  much  of 
the  kinetic  energy  in  a  damped  system  could  be  dissipated  before  being  transmitted.  Depending 
on  the  parameters  of  the  system,  the  response  may  exhibit  underdamped  beats  as  shown  in  Fig¬ 
ure  2,  critically  damped  beats,  or  overdamped  beats.  In  order  to  demonstrate  the  physical  and 
analytical  basis  of  the  method  we  consider  the  special  case  of  a  damped  tuned  system  with 
y  >  (p  -  B„)2.  Other  cases  can  be  handled  in  exactly  the  same  manner. 

When  |=0  and  y  >  (p-B„)2  the  contribution  from  the  tuning  poles  in  eq?  (21) 
becomes 

w(r)  =  -  7 - - r —  f  m,(t)  [r  T)  2C0Sa,(^ — r)  sinTjfr— t )] z/t  (27) 

[y-(p-Bn)2]h  J0 


where 

V  =  [y-  (p-  Bn)2]'l:<ol2 
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We  note  that  the  term  in  brackets  in  the  integral  represents  a  damped  beat  motion  as  shown  in 
Figure  2  with  a  damping  factor  (/3  +  B„)/2  and  a  beat  period  T  given  by  2n/r).  When  the  term 
sinij(r-r)  is  expanded 


ii(t)  =  - 


+ 


C> 


[y-Q3-fl„)2l* 


-  cos(r}r-0) 
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where 


9  =  tan  1 
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We  consider  r\t\  =  2irt\lT  «  1,  where  t|  is  the  duration  of  ground  motion  and  7"  is  the 
beat  period  of  the  system.  Then  the  first  integral  in  Eq.  (28)  can  be  approximated  by 

J  Ug\j)  e  ”  cosojW-tMt 

o 


and  the  second  neglected  since  sin*}/  will  be  bounded  by  r}<\  «  1  and  uK  -  0  for  t  >  tt. 
Thus,  we  take 


’u(t)  =  - 


C"(>>  sin-rjt 

T^Os-fl,,)2]" 


J*  Ug(r)  e  r)/2 cosoj(i~t)(Jt 


(30) 


When  the  parameters  y'\  /3,  and  Bn  are  small,  this  result  may  be  interpreted  in  the  fol¬ 
lowing  way:  for  t  >  th  the  above  expression  can  be  written  in  the  form 

C"cj2s\ni)t  -[fj+B 


ii(t)  =  - 


where 


with 


2tj 

R  =  (. 4?+A})V: 


R  cos  ( wt-il >) 


.  t  ■■  t  \  +U}+B„)wl/2  , 

A  |  =  J  ug(t)  e  '  cos wt  dt 

o 

.  t  ■■  t  \  +lP+BJ»l/2  .  ,  , 

A 2  —  J  Ugil)  e  *  sin mt  dt 
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and 

0  =  tan-1  (Aii  A\) 


The  response  indicated  by  the  above  is  illustrated  in  Fig.  5.  In  the  above,  the  terms  R  and  < // 
are  constants  independent  of  t  for  t  >  r (,  and  R  cos(a>/-</<)  is  a  rapidly  varying  function  of 
time.  The  term 


or  -(/H#  )ui/2  . 
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is  a  slowly  varying  envelope  curve  whose  maximum  value  must  be  determined.  The  maximum 
value  of  this  envelope  curve  is  attained  at  time  /*,  expressed  by 


tan  i)t*  = 


(fi+B„U 


(31) 


The  value  of  t*  is  thus 


t*  =  aretan[2-r)/<D  (($+B„)]/ri 


(32) 


For  lightly  damped  systems  and  light  equipment,  in  general  t*  »  tx.  The  values  of  sinrjt  and 
expt-(j8+^fl)<D//2]  when  the  maximum  of  the  envelope  is  achieved  are 


simjt*  = 


[t)2+(/3+5„)2oj2/4]'/! 


-(/J+fl„)wfV2 

e  " 


(33) 


where 

k  =  (arctanO/i 
{  =  (y-(fi-B„)2]'h/(p+B„) 

It  follows  that 

I  &  I  max  = 

r* 

(p+B„)a> r,f/2.  f  ■■  ,  v  -(B+BJojU'- t)/2  .  ,  ,  ,  .  ,, 

e  | J  ug(r)e  "  cosoj(t*-r)</r|  (36) 

o 

In  order  that  this  estimate  of  peak  acceleration  be  useful  for  design  purposes,  it  is  neces¬ 
sary  that  the  second  factor  in  braces  be  interpreted  in  terms  of  a  ground  response  spectrum.  To 
this  end,  we  recognize  that  the  integral  is,  to  the  order  of  P+B„,  the  relative  velocity  response 
history  evaluated  at  time  t*  of  a  lightly  damped  single-degree-of-freedom  oscillator  of  frequency 
to  and  damping  factor  (/3+B„)/2  subjected  to  the  ground  acceleration  ug(r).  At  some  time  t 
during  the  ground  motion  or  shortly  after  it  ceases  (so  that  i  «  r*),  the  absolute  value  of  the 
relative  velocity  will  attain  its  global  maximum,  denoted  as  |v(r)|.  The  relative  velocity 
response  at  t*  denoted  as  v(r*),  can  be  thought  of  as  that  which  would  occur  in  a  single- 
degree-of-freedom  system  (subjected  to  the  ground  acceleration  ug(t))  as  a  consequence  of  free 
vibration  beginning  at  time  /(>/,)  when  the  absolute  value  of  the  relative  velocity  of  the  oscil¬ 
lator  attains  its  first  JocaJ  maximum,  |  v  ( r )  | ,  after  the  end  of  the  ground  motion.  This  instant  of 
time  t  is  equal  to  7  if  t  occurs  after  the  end  of  the  ground  motion;  otherwise,  t  >  l  In  any 
event,  t  «  t*.  Thus,  we  can  write 

|  v  (/*)  |  =  |  v  (/)  |e~</3+B">w<,"-,)/2|cosci>(r*— /)  |  (37) 

where  |v(r)|  <  |v(/)|.  It  then  follows  that 
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since  t/t*  «  1.  From  this  we  obtain  the  approximate  result 


| v (F) |  =  e(f3+B")m,'l2\  J  ug(r)e 


-(B+B„)u(r‘-T)/2 


cosw(f*-r)rfT| 


(39) 


We  recognize,  however,  that  to  the  order  of  /3  and  B„,  |  v(F)  |  is  very  nearly  the  pseudo- velocity 
response  spectrum  S^(<o,  (fi+B„)/ 2)  for  a  lightly  damped  single-degree-of-freedom  oscillator  of 
frequency  o»  and  damping  factor  (fi+B„)/2  subjected  to  the  ground  acceleration  iitU).  Thus,  an 
estimate  of  the  maximum  equipment  acceleration  is 


\Crn\m\sinr)t*\e~*  P+Bn 
[y-(fi-BnW>  Syi"’  2 


With  the  value  of  sin-rjt*  from  Eq.  (33),  we  obtain  the  final  estimate  as 

I" I™*  “  (y+4/3 B„y/! 


(40) 


Recalling  that 


0)Sy  8=1  SA 


“  w2Sd 


this  estimate  can  be  written  in  the  alternative  form 


|  it  ( t )  |  max 
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SA(ai, 
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(41) 


Estimates  for  other  cases  of  damped  tuned  systems,  e.g.,  where  y  <  (/3  -  Bn)2,  y  -  (/3  -  B„)2 
and  y  -  B„ 2,  /3  -  0,  were  developed  in  [111;  precisely  the  same  result  was  obtained. 

For  a  slightly  nontuned  system,  the  dominant  contribution  to  the  response  is  given  by 
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where  X  and  fi  are  defined  in  terms  of  y ,  /3,  and  B„  in  Eqs.  (16)  and  (17). 

The  reduction  of  this  general  form  to  the  required  amplified  response  spectrum  form  is 
given  in  detail  in  reference  1 1  and  takes  the  form 
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The  contribution  to  the  equipment  acceleration  from  the  other  modes  for  which  no 
interaction  need  be  considered  is  quite  standard.  Using  eqs.  (8)  and  (9)  in  conjunction  with  the 
square  root  of  the  sum  of  squares  procedure  the  estimate  of  the  maximum  acceleration  in  the 
early  peak,  which  will  occur  during  the  ground  motion,  is  given  by 
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The  peak  value  given  by  eq.  (46)  occurs  early  after  only  a  few  cycles  of  equipment  motion 
and  that  given  by  eq.  (43)  occurs  late  and  thus  in  principle  they  should  not  be  taken  together. 
However,  we  have  recently  completed  extensive  numerical  and  physical  experimentation  the 
results  of  which  appear  to  suggest  ways  in  which  the  results  can  be  conveniently  applied  to 
equipment  mounting  design.  The  physical  experiments  were  carried  out  at  the  Earthquake 
Simulator  Laboratory  of  the  Earthquake  Engineering  Research  Center,  University  of  California, 
Berkeley,  and  involved  the  use  of  a  three-story  steel  frame  model  structure  (Figure  3)  around 
one-third  full  scale.  The  response  of  the  steel  frame  model  is  roughly  that  of  a  shear  building 
with  a  natural  frequencies  at  around  2  Hz,  8  Hz,  and  15  Hz.  The  total  weight  of  the  model  *was 
39,500  lbs  and  its  height  was  about  20  feet.  Three  single-degree-of-freedom  oscillators  (Figure 
4)  were  attached  to  the  concrete  blocks  at  the  second  and  third  floors  to  simulate  equipment  in 
a  primary  structure.  The  test  structure  was  instrumented  to  measure  displacements  and 
accelerations  at  each  floor  and.  acceleration  of  the  oscillators.  The  mechanical  oscillators,  each 
weighing  about  20  lbs,  were  constructed  to  correspond  to  the  first  three  natural  frequencies  of 
the  model  structure. 

The  first  three  natural  frequencies  of  the  frame  were  determined  accurately  and  the  oscil¬ 
lators  were  tuned  to  the  structural  frequencies.  They  were  then  bolted  to  the  floors  of  the 
model  structure.  A  variety  of  ground  shock  table  motions  was  applied  to  the  model,  some  of 


which  were  short  duration  approximate  square  waves  and  some  of  which  were  longer  duration 
and  more  erratic.  Figure  5  shows  the  time  history  of  acceleration  in  a  typical  run  for  a  short 
duration  single  displacement  pulse.  This  produced  peak  table  acceleration  pulses  around 
±0.16g,  separated  by  about  1  sec.  The  peak  responses  in  the  tuned  oscillators  were  2.0g,  3.5g, 
and  0.6g  in  the  first,  second,  and  third  mode  oscillators,  respectively.  On  Figure  5,  the  top 
trace  shows  the  input  table  acceleration,  the  second  the  response  of  the  third  floor  of  the  struc¬ 
ture,  and  the  next  three  the  accelerations  of  the  three  tuned  oscillators.  The  very  large 
magnification  of  acceleration  experienced  by  the  oscillators  is  immediately  obvious.  The  beat 
phenomenon  is  clear.  The  peak  acceleration  in  the  first  two  oscillators  is  achieved  considerably 
after  the  peak  acceleration  in  the  input.  These  two  oscillators  are  obviously  responding  at  the 
coupled  frequencies  governed  by  equipment-structure  interaction. 

The  response  of  the  third  oscillator  is  quite  different;  a  highly  irregular  pattern  appears 
and  local  maxima  occur  during  as  well  as  after  the  excitation.  Fourier  transforms  of  the 
acceleration  time  histories  of  the  three  oscillators  have  been  taken  to  clarify  this  response  (Fig¬ 
ure  6).  It  is  clear  from  these  transforms  that  the  first  mode  oscillator  responded  only  at  the 
first  mode  frequency,  around  2  Hz;  the  second  mode  oscillator  responded  predominantly  at  the 
second  mode  frequency,  around  8  Hz,  with  a  small  contribution  from  the  first  mode.  The 
dominant  contribution  to  the  response  of  the  third  mode  oscillator  was  from  frequencies  around 
the  third  mode  frequency  of  15  Hz,  but  significant  contributions  also  appeared  from  the  lower 
modes.  Thus,  the  pure  damped  beat  response  predicted  by  the  theory  was  complicated  by  con¬ 
tributions  from  the  lower  modes. 

These  observations  are  supported  by  the  numerical  experiments  and  have  implications  for 
the  application  of  the  theoretical  results  to  design.  First,  only  the  structural  modes  with  fre¬ 
quencies  up  to  and  around  that  of  the  equipment  under  consideration  need  be  considered  and, 
second,  the  late  peak  acceleration  given  by  eq.  (43)  should  be  summed  with  the  early  peak 
acceleration  given  by  eq.  (46)  by  an  appropriate  summation  rule.  For  instance,  if  the  conven¬ 
tional  square  root  of  the  sum  of  the  squares  procedure  is  used,  the  estimate  is 
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If  the  late  peak  occurs  well  after  the  peak  of  the  excitation  (as  in  the  first  mode  oscillator)  it 
will  be  the  dominant  term  and  the  result  will  be  nearly  the  same  as  if  it  alone  were  considered. 
On  the  other  hand,  if  the  late  peak  occurs  during  the  excitation  (as  in  the  third  mode  oscillator) 
then  it  should  be  superposed  with  the  other  modal  contributions  which  may  occur  around  the 
same  time  and  this  will  be  effected  by  the  formula. 


NUMERICAL  EXPERIMENTS 

In  addition  to  the  physical  experiments,  we  have  performed  numerical  experiments  using 
a  standard  structural  analysis  program,  TABS  [12],  to  compute  by  a  variety  of  methods  the 
response  of  a  light  appendage  in  a  building.  The  structure  used  in  these  numerical  experiments 
was  a  ten-story  reinforced  concrete  frame  building  described  and  used  in  reference  13;  the 
appendage  was  a  single-degree-of-freedom  oscillator  attached  to  the  top  floor  and  had  a  mass 
which  gave  a  mass  ratio  of  0.001  compared  to  the  modal  mass  of  the  first  mode  of  the  building. 
The  response  of  the  appendage  was  calculated  using  TABS  to  evaluate  the  eleven-degree-of- 
freedom  system  which  consisted  of  the  structure  and  the  appendage.  The  program  has  the 
capability  of  computing  natural  frequencies  and  mode  shapes  and  then  can  either  compute  a 


ime  history  in  each  mode  (by  a  method  which  is  exact  for  piecewise-linear  acceleration 
records)  and  sums  these  for  the  resultant  response,  or  can  determine  the  maximum  response  in 
each  mode  by  a  spectrum  method  and  estimate  the  maximum  response  by  the  square  root  of 
the  sum  of  the  squares  procedure. 

The  response  of  the  appendage  as  a  function  of  appendage  frequency  (maintaining  the 
same  mass  ratio)  was  calculated  by  both  methods  for  a  variety  of  ground  motions.  Two  cases  of 
damping  were  considered:  undamped  and  2%  of  critical  damping  in  each  mode.  Typical  results 
for  a  very  short  duration  ground  motion  for  both  methods  of  computation  are  shown  in  Figure 
7  for  the  undamped  case  and  in  Figure  8  for  the  damped  case.  The  TABS  time-history  results, 
which  should  be  exact,  are  shown  by  solid  dots  and  the  TABS  spectrum  (SRSS)  results  are 
shown  by  triangles.  The  tuned  system  estimate  based  on  eq.  (47)  is  shown  as  a  solid  curve  and 
the  nontuned  system  estimate  based  on  eq.  (24)  is  indicated  by  squares  on  the  same  figures. 


It  is  clear  from  the  diagrams  that  for  both  damped  and  undamped  cases  the  tuned  esti¬ 
mate  compares  well  with  the  nontuned  estimate  in  the  regions  of  gross  nontuning.  (This  is 
probably  due  to  the  even  spacing  of  the  natural  frequencies  of  the  structure  in  this  case;  thus, 
the  tuned  estimate  could  be  used  for  structures  with  relatively  evenly  spaced  natural  frequen¬ 
cies  except  possibly  for  frequencies  well  below  the  first  or  much  above  the  highest.)  Further, 
the  tuned  estimate  compares  well  with  the  time-history  calculation  over  the  entire  range  of 
appendage  frequencies  considered,  which  encompasses  the  first  three  modes  of  the  structure. 
This  is  surprising  in  view  of  the  fact  that  the  spectra  used  in  these  numerical  experiments  are 
jagged  spectra  and  not  the  smoothed  design  spectra  for  which  the  approach  was  developed.  A 
further  result  is  that  the  TABS  spectrum  (SRSS)  estimate  is  very  poor  near  tuning  for  all 
modes  considered,  but  is  very  good  for  grossly  nontuned  cases.  This  is  not  surprising  since  it  is 
well  known  that  SRSS  is  not  accurate  for  closely  spaced  modes  which  naturally  occur  when  a 
light  appendage  or  equipment  has  a  frequency  near  one  of  the  structure’s  frequencies.  In  this 
case,  the  standard  approach  is  to  use  the  NRL  method,  which  for  a  two-degree-of-freedom  sys¬ 
tem  is  equivalent  to  the  absolute  sum  method.  From  Figure  7,  the  TABS  time-history  result 
for  the  undamped  case  and  when  the  appendage  is  tuned  to  the  first  mode  is  78.2.  The  TABS 
spectrum  (SRSS)  result  is  55.1  and  if  the  absolute  sum  method  is  used  the  result  is  77.9.  The 
estimate  given  by  the  formula  (47)  for  the  tuned  case  is  78.1.  However,  the  TABS  spectrum 
(SRSS)  result  greatly  overestimates  the  response  in  the  damped  case  as  shown  in  Figure  8.  The 
TABS  time-history  computation  is  20.2;  the  TABS  spectrum  (SRSS)  result  is  50.2.  (If  the 
absolute  sum  method  is  used,  the  result  is  71.0.)  The  tuned  result  from  eq.  (47)  is  19.4. 


The  physical  explanation  for  this  overestimation  by  the  SRSS  and  absolute  sum  methods 
in  the  damped  case  is  that  the  beat  phenomenon  involves  an  energy  transfer  between  the  two 
elements  of  the  system.  In  the  case  of  an  undamped,  perfectly  tuned  system,  the  energy 
transfer  is  complete,  i.e.,  maximum  response  in  one  element  is  accompanied  by  zero  response 
in  the  other.  This  energy  transfer  takes  time;  the  beat  period  that  controls  the  energy  transfer 
is  inversely  proportional  to  the  difference  between  the  two  closely  spaced  frequencies.  The 
phenomenon  can  be  interpreted  geometrically  in  terms  of  the  eigenvectors  of  the  modes.  The 
components  of  these  modes  can  be  thought  of  as  vectors  in  a  generalized  state  space.  The 
components  that  represent  the  appendage  are  initially  1 80  degrees  out  of  phase  and  those  of  the 
structure  are  in  phase.  As  the  motion  continues,  the  equipment  components  rotate  and  eventu¬ 
ally  align,  while  the  structural  components  become  out  of  phase.  As  the  vectors  rotate,  the 
resultant  for  the  equipment  increases  and  if  the  system  is  undamped,  attains  a  maximum  when 
they  are  aligned.  In  the  damped  case,  the  peak  values  of  the  vectors  will  first  increase  and  then 
diminish  as  they  rotate  and  the  resultant  will  achieve  its  maximum  before  they  line  up.  Thus, 
it  is  important  in  the  damped  case  to  determine  at  which  point  the  maximum  of  the  resultant 
occurs.  In  terms  of  this  analogy,  the  square  root  of  the  sum  of  the  squares  procedure  assumes 
that  the  peak  values  of  both  of  the  vectors  are  attained  simultaneously  and  occur  when  they  are 
at  90  degrees.  The  absolute  sum  method  assumes  that  this  occurs  when  they  are  lined  up.  The 
method  described  here  evaluates  the  position  at  which  the  resultant  actually  attains  its  max¬ 
imum. 
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The  overestimation  at  tuning  of  the  SRSS  or  absolute  sum  methods  for  cases  with  damp¬ 
ing  can  be  quantified  by  considering  only  the  two-degree-of-freedom  system  equivalent  to  the 
tuned  appendage-structure  system.  The  three  estimates  of  peak  acceleration  based  on  the  SRSS 
rule,  the  absolute  method,  and  our  formula,  eq.  (47),  respectively,  are  given  by 

M  =  ^  sA<o,p) 

where  for  simplicity  B„  is  taken  to  be  the  same  as  p.  We  note  immediately  that  if  p  — *0,  the 
square  root  of  the  sum  of  the  squares  estimate  will  be  low  by  a  factor  of  V2;  the  absolute  sum 
estimate  and  the  new  estimate  are  in  this  case  identical.  However,  the  estimates  differ  for 
nonzero  damping.  We  introduce  an  overestimation  ratio 

R  =  f-7222-  =  (1  +  4/3  2/y)  V  (48) 

y  new  ^  ' 

!•*  I  max 


Since  k  is  given  by 

k  =  arctan 

the  overestimation  ratio  can  be  expressed  in  terms  of  the  single  parameter  yVl/2p,  in  terms  of 
which  a  plot  of  eq.  (48)  is  given  in  Figure  9.  The  result  can  also  be  shown  for  several  values  of 
yasa  function  of  /3;  the  resulting  curves  are  shown  in  Figure  10.  Clearly,  for  all  values  of  y, 
the  overestimation  parameter  and  thus  the  conservatism  of  the  conventional  method  steadily 
increase  with  p.  If  we  fix  y  and  decrease  P  such  that 

yV:/2p  — *  00 


yv,/2p 

yh!2p 


we  find  that 


k  — > 


TT  2P 

2  y'/! 


and 

*  -  „  +  SLu ■ 

v 

which  indicates  that  R  >  1  for  all  nonzero  p. 

On  the  other  hand,  if  we  fix  p  and  decrease  y  such  that  yl/V2/3— 0,  we  find  that  for  light 
appendage  cases,  which  correspond  to  very  closely  spaced  modes,  the  response  can  be  greatly 
overestimated. 

It  is  interesting  to  compare  these  results  to  those  from  the  numerical  experiments  on  the 
multidegree-of-freedom  structure.  The  absolute  sum  method  estimate  for  the  2%  damped  case 
(as  shown  in  Figure  8)  is  71.0.  The  new  estimate,  from  eq.  (47),  is  19.4,  and  the  time-history 
result  is  20.2.  The  overestimate  of  the  absolute  sum  method  is  thus  71.0/20.2  —  3.52;  for  y  “ 
0.001  and  p  —  0.02,  the  value  of  R  is  3.76.  This  indicates  the  excellent  applicability  of  the 
simple  analysis  based  on  the  one-degree-of-freedom  structure  to  the  multidegree-of-freedom 
structure.  The  reason  for  this  is  the  dominance  of  the  two  closely  spaced  modes  in  the  compo¬ 
site  system  when  tuned. 
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OPTIMUM  DAMPING 

The  solution  for  the  nearly  tuned  two-degree-of-freedom  system  given  in  an  earlier  sec¬ 
tion  depends  on  a  considerable  number  of  parameters  and  it  is  difficult  to  interpret  the 
influence  of  individual  terms.  The  damping  can  obviously  dominate  the  solution  if 
4/8 B  »  y  +  £2,  but  it  is  also  clear  that  the  difference  between  damping  factors  in  the  structure 
and  the  equipment  can  play  a  significant  role  in  the  term  £  and  thus  in  the  amplification  factor. 
Thus,  it  is  appropriate  to  ask  the  question:  for  what  values  of  p  and  B,  given  fixed  values  of  y, 
£2,  and  p  +  B,  is  the  amplification  A  minimized?  To  this  end,  we  denote 

y  +  £2  =  a2 


P  +  B  =  c 

and  take  a  and  c  as  fixed  positive  numbers. 

We  also  note  that  in  the  case  where  y  +  £2  >  (/3 -B)2,  the  denominator  (y  +  £2  +  4/3 B)Vl 
can  be  written  in  the  form  c(l  +  £)'/!.  The  amplification  function  A  becomes 

e~x(o 

A  =  c  (1  +  £2)'/: 

Considered  as  a  function  of  £,  A  has  extremal  values  when 


(l  +  £2)'/!  di 


(l+£2): 


dK  1  ,  ..  .  1  1 


and  it  follows  that 


1  di  £(l+£2)'/!  |  i  1-K2  1+£2J  £(1+£2)i/!  '  £ 

Thus,  dA/dl,  — *  0  when  £—*<».  When  £  — •  0,  we  have 

cd±  (1  _  1  {2  _  d 

dC  £  3  s 

There  are  no  other  zeroes  of  dA/d£.  We  note  that  as  £  — »  0,  A  (£)  has  the  form 

/!(£)  =  — (1  -  £) 

c  b 

Thus,  the  function  A  (£)  is  a  monotonically  decreasing  function  of  £.  To  minimize  A,  then,  it 
is  necessary  to  maximize  £  and  this  is  achieved  obviously  when  /3  —  B.  The  conclusion  is 
therefore  that  when 

y  +  £2  >  (p-  B)2 

the  optimal  damping  for  fixed,  p  +  B  is  p  -  B.  The  value  of  A  for  this  case  is  given  by 

A  — - — - 

mm  (y +  £2+  c2)/i 

with  k  -  arctan  £/£  and  £  -  (y-t -f2)'/:/c.  We  note  that  this  value  is  always  less  than  e~K/c.  If  c 
is  sufficiently  large,  it  is  also  possible  to  have 
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y  +  f  <  (p  -  B)2 


and  in  this  case  the  amplification  factor  takes  the  form 

j  « - — - 

(y  +  £2  +  4/35),/: 

with  k  —  arctanh  £/£,  where  now 

Up-B)2- y-Q)‘/; 
4  “  £ 


In  this  case  the  denominator 


(y  +  S2  +  ApB)'h  =  c(l-£2)'/! 


Thus, 


C  A  (£> 


„-*•{) 


{\-l2)'h 

Using  the  well-known  expression  for  arctanh  in  terms  of  In,  we  have 


c 


ii-O 

(l-C)  { 


(l+{) 

(l+o  { 


and 


Thus 


In  (cA )  -  -^-{(1  -£)ln(l  -0  “  <1  +  0  In  (I  +0) 
2d{^A)  =  -L  (_in(1  _  f)  +  In (1  +  0  -  2£) 


Thus,  dA/dt,  — *  0  when  In  1(1+0/0-01  '=  2£  if  A  &  This  has  only  one  solution,  namely 
£  -0. 

The  function  cA  (0  has  two  parts,  e~K  and  l/(l-£2)'/!.  When  £  — 1 ►  0, 
e~*  — *  e_1(l  -  l/3£2),  but  1/(1— £2)'/!  —  1  + 1/6£2.  Thus,  — 1 >l  +  l/6£2  for  small  £  and 

approaches  a  limiting  value  as  £  — ♦  1;  the  limit  of  A  as  £  — *  1  is  given  by 

lim cA  _  lim 

£-1  "£-l 


JdL 

(1-0  1 


1+L 

(1  +£)  { 


Since  the  limit  of  cA  as  £  —  0  is  e~'  and  dA/d£  =  0  only  at  £  =  0,  the  function  increases 
monotonically  from  e~'  at  £  —  0,  to  1/2  at  £  =  1.  The  minimum  value  of  A  is  achieved  here  by 
making  £  as  small  as  possible.  If  P  and  B  are  such  that  (y  +  £2)  can  be  equal  to  (p-  B)2,  then 
the  minimum  is  at  £  -  0  and  has  the  value  e~l/c,  and  P  and  fiare  given  by 

p  =  |  (c  ±  (y  +  m 

B  -  |  (c  +  (y  +  €2)Vl) 


This  requires  that  c  >  (y  +  £2)'/j  in  order  to  have  realistic  damping,  but  in  any  case  the  analysis 
of  the  previous  case,  y +  £2  >  (p-  B)2,  gives  a  minimum  for  A  which  is  always  less  than  e~'/c. 
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It  follows  then  that  in  all  cases  the  optimum  choice  of  and  B  is  /3  —  B. 

CONCLUDING  REMARKS 

The  advantages  of  the  estimation  technique  developed  in  this  paper  are  its  simplicity  and 
adaptability  for  practical  application.  A  great  deal  of  computational  effort  is  avoided  since 
time-history  analyses  need  not  be  performed.  The  equipment  and  structure  need  not  be 
analyzed  as  an  N+l  -degree-of-f reedom  system  either  by  modal  or  matrix-time-marching 
methods,  and  errors  in  estimates  of  peak  response  due  to  the  possible  unreliability  of  numerical 
time  integration  schemes,  or  to  uncertainty  as  to  the  appropriate  procedure  for  summing  the 
contributions  of  the  two  closely  spaced  modes,  are  thereby  avoided.  For  tuned  and  nearly 
tuned  systems,  the  method  accounts  for  the  important  effect  of  equipment-structure  interac¬ 
tion,  which  effect  is  completely  neglected  in  the  floor  spectrum  method.  The  method  advanced 
here  does  not  require  that  new  information  be  generated.  Data  available  from  the  building 
design  alone,  the  equipment  design  alone,  and  the  response  spectra  are  used.  The  estimates  of 
peak  response  have  been  obtained  by  rational  analysis  and  are  easily  evaluated  and  conveniently 
used  during  the  design  process. 

There  are  a  number  of  obvious  extensions  of  the  analysis  on  which  these  simple  estimates 
are  based  which  could  be  of  significance  to  DNA  problems.  The  method  which  in  effect  intro¬ 
duces  a  transfer  function  between  the  induced  ground  motion  and  the  response  of  the  equip¬ 
ment  avoids  the  explicit  assumption  of  the  support  motion.  For  example,  equipment  such  as  a 
missle  in  a  protective  structure  and  suspended  in  this  structure  at  several  points  could  be 
treated  by  this  approach.  In  the  standard  approach  to  multiple  support  response  spectrum 
analysis  of  such  a  system  the  motion  of  each  support  would  be  calculated  as  if  the  suspended 
system  were  not  present  and  this  support  motion  applied  to  the  suspended  item.  To  use  a 
response  spectra  method  for  such  a  system  it  would  be  necessary  either  to  envelop  the  response 
spectra  of  the  different  support  points  or  else  to  assume  that  each  support  point  moves  indepen¬ 
dently  of  the  others  and  superpose  the  responses  by  a  summation  rule.  This  must  provide  a 
highly  conservative  prediction  since  it  ignores  the  fact  that  all  support  motions  are  produced  by 
a  single  input.  It  further  cannot  take  into  account  closely  spaced  modes  in  the  suspended  item 
and  the  protective  structure  and  this  will  provide  even  more  conservatism  as  outlined  in  the 
paper.  A  straightforward  extension  of  the  present  approach  will  provide  response  spectra  esti¬ 
mates  of  the  motion  of  the  suspended  item  in  terms  of  the  ground  motion  utilizing  the  transfer 
functions  for  the  support  points,  but  not  explicitly  calculating  their  motion.  Work  on  this  is 
continuing. 

A  further  extension  of  the  analysis  which  is  possible  and  could  apply  to  DNA  problems  is 
to  very  long  duration  inputs  which  are  known  only  in  a  statistical  sense.  Here,  it  is  not  possible 
to  exploit  the  beat  phenomenon  and  the  consequent  late  peak  to  obtain  the  response  spectrum 
result,  but  the  transfer  function  result  used  as  a  point  of  departure  still  holds.  Using  this  and 
the  other  assumptions  of  the  analysis  it  will  be  possible  to  obtain  statistical  estimates  of  peak 
equipment  response  in  terms  of  amplification  factors  and  the  statistics  of  the  input  ground 
motion.  This  could  be  a  particularly  valuable  approach  to  the  problem  of  a  transporter  vehicle 
carrying  a  sensitive  item  of  equipment  over  a  road  or  track  with  a  statistically  specified  surface 
roughness. 
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Maximum  Equipment  Displacement,  in  Feet 


—  Tuned  Estimate 
°  Untuned  Estimate 


•  TABS -Time  History 


a>|  u>2 

ca,  Equipment  Frequency,  in  Radians/Sec. 


Appendage  Respone — Short  Duration  Shock — 2%  Damping 


26 


REFERENCES 


[1]  J.  M.  Kelly  and  J.  L.  Sackman,  "Response  Spectra  Design  Methods  for  Tuned 
Equipment-Structure  Systems,"  Journal  of  Sound  and  Vibration,  Vol.  59,  No.  2,  pp.  171- 
179  (1978). 

[2]  J.  M.  Kelly  and  J.  L.  Sackman,  "Shock  Spectra  Design  Methods  for  Equipment-Structure 
Systems,"  to  appear  in  Bulletin  of  Shock  and  Vibration  (1979). 

[3]  N-C.  Tsai,  "Spectrum-Compatible  Motions  for  Design  Purpose,"  Journal  of  the  Engineering 
Mechanics  Division,  Proceedings  ASCE,  Vol.  98,  EM2,  pp.  345-356  (1972). 

[4]  J.  M.  Biggs  and  J.  M.  Roesset,  "Seismic  Analysis  of  Equipment  Mounted  on  a  Massive 
Structure,"  in  Seismic  Design  for  Nuclear  Power  Plants,  edited  by  R.  J.  Hansen,  MIT  Press, 
Cambridge,  Massachusetts  (1970). 

[5]  K.  K.  Kapur  and  L.  C.  Shao,  "Generation  of  Seismic  Floor  Response  Spectra  for  Equip¬ 
ment  Design,"  Proceedings,  Structural  Design  of  Nuclear  Power  Plant  Facilities,  Chicago, 
Dec.  17-18,  1973,  American  Society  of  Civil  Engineers,  New  York,  N.  Y.  (1973). 

[6]  M.  Amin,  et  ai,  "Earthquake  Response  of  Multiply  Connected  Light  Secondary  Systems 
by  Spectrum  Methods,"  Symposium  on  Seismic  Analysis  of  Pressure  Vessel  and  Piping  Com¬ 
ponents,  ASME,  New  York,  pp.  103-129  (1971). 

[7]  N.  M.  Newmark,  "Earthquake  Response  Analysis  of  Reactor  Structures,"  Nuclear  Engineer¬ 
ing  and  Design,  Vol.  20,  No.  2,  pp.  303-322  (1972). 

[8]  A.  H.  Hadjian,  "Some  Problems  with  the  Calculation  of  Seismic  Forces  on  Equipment," 
Proceedings,  Structural  Design  of  Nuclear  Plant  Facilities,  Chicago,  Dec.  17-18,  1973, 
American  Society  of  Civil  Engineers,  New  York,  N.  Y.  (1973). 

[9]  J.  Penzien  and  A.  K.  Chopra,  "Earthquake  Response  of  an  Appendage  in  Multi-Story 
Building,"  Proceedings,  Third  World  Conference  on  Earthquake  Engineering,  Vol.  II,  New 
Zealand  (1965). 

[10]  K.  A.  Peters,  D.  Schmitz,  and  U.  Wagner,  "Determination  of  Floor  Response  Spectra  on 
the  Basis  of  the  Response  Spectrum  Method,"  Nuclear  Engineering  and  Design,  Vol.  44.  pp. 
255-262  (1977). 

[11]  J.  L.  Sackman  and  J.  M.  Kelly,  "Rational  Design  Methods  for  Light  Equipment  in  Struc¬ 
tures  Subjected  to  Ground  Motion,"  Report  No.  UCB/EERC-78/19,  Earthquake  Engineer¬ 
ing  Research  Center,  University  of  California,  Berkeley  (1978). 

[12]  E.  L.  Wilson  and  H.  Dovey,  "Three-Dimensional  Analysis  of  Building  Systems,"  Report 
No.  EERC  72-8,  Earthquake  Engineering  Research  Center,  University  of  California, 
Berkeley  (1972). 

[13]  S.  W.  Zagajeski  and  V.  V.  Bertero,  "Computer-Aided  Optimum  Design  of  Ductile  Rein¬ 
forced  Concrete  Moment-Resisting  Frames,"  Report  No.  UCBIEERC-77116,  Earthquake 
Engineering  Research  Center,  University  of  California,  Berkeley  (1977). 


29 


DISTRIBUTION  LIST 


DEPARTMENT  OF  DEFENSE 

Defense  Intelligence  Agency 
ATTN:  RDS-3A 
ATTN:  DB-4C2,  B.  Morris 

Defense  Nuclear  Agency 
ATTN:  SPAS 
ATTN:  STSP 
ATTN:  SPSS 
4  cy  ATTN:  TITL 

Defense  Technical  Information  Center 
12  cy  ATTN:  DD 

Field  Command 
Defense  Nuclear  Agency 
ATTN:  FCPR 
ATTN:  FCT 
ATTN:  FCTMOF 

Joint  Strat.  Tgt.  Planning  Staff 
ATTN:  DOXT 
ATTN:  XPFS 
ATTN:  JLTW-2 
ATTN:  NRI-STINFO  Library 

NATO  School  (SHAPE) 

ATTN:  U.S.  Documents  Officer 

DEPARTMENT  OF  THE  ARMY 

BMD  Advanced  Technology  Center 
Department  of  the  Army 
ATTN:  1CRDABH-X 
ATTN:  ATC-T 

BMD  Systems  Command 
Department  of  the  Army 

ATTN:  BMDSC-H,  N.  Hurst 

Chief  of  Engineers 
Department  of  the  Army 
ATTN:  DAEN-MCE-D 
ATTN:  OAEN-RDM 

Construction  Engineering  Rsch.  Lab. 

Department  of  the  Army 
ATTN:  CERL-SOI-L 

Engineer  Studies  Center 
Department  of  the  Army 

ATTN:  DAEN-FES,  LTC  Hatch 

Harry  Diamond  Laboratories 
Department  of  the  Army 
ATTN:  DELHD-N-P 
ATTN:  DELHD-I-TL 

U.S.  Army  Armament  Material  Readiness  Command 
ATTN:  MA,  Library 

U.S.  Army  Ballistic  Research  Labs. 

ATTN:  DRDAR-BLV 

ATTN:  DROAR-BLE,  J.  Keefer 


DEPARTMENT  OF  THE  ARMY  (Continued) 

U.S.  Army  Communications  Command 

ATTN:  Technical  Reference  Division 

U.S.  Army  Engineer  Center 
ATTN:  ATZA 

U.S.  Army  Engineer  Div.,  Huntsville 
ATTN:  HNDEL  SR 

U.S.  Army  Engineer  Div.,  Ohio  River 
ATTN:  ORDAS-L 

U.S.  Army  Engineer  School 
ATTN:  ATZA-DTE-ADM 
ATTN:  ATZA-CDC 

U.S.  Army  Engr.  Waterways  Exper.  Station 
ATTN:  WESSS,  J.  Ballard 
ATTN:  Library 
ATTN:  WESSD,  G.  Jackson 
ATTN:  WESSE,  L.  Ingram 
ATTN:  J.  Strange 
ATTN:  WESSA,  W.  Flathau 

U.S.  Army  Foreign  Science  A  Tech.  Ctr. 

ATTN:  DRXST-SD 

U.S.  Army  Mat.  Cmd.  Proj.  Mngr.  for  Nuc.  Munitions 
ATTN:  DRCPM-NUC 

U.S.  Army  Material  &  Mechanics  Rsch.  Ctr. 

ATTN:  Technical  Library 

U.S.  Army  Materiel  Dev.  &  Readiness  Cmd. 

ATTN:  DRXAM-TL 

U.S.  Army  Missile  R&D  Command 
ATTN:  DRDMI-XS 
ATTN:  RSIC 

U.S.  Army  Mobility  Equip.  R&D  Cmd. 

ATTN:  DRDME-WC 

U.S.  Army  Nuclear  &  Chemical  Agency 
ATTN:  Library 

U.S.  Army  War  College 
ATTN:  Library 

U.S.  Military  Academy 

ATTN:  R.  La  Frenz 

DEPARTMENT  OF  THE  NAVY 

Marine  Corps 

Department  of  the  Navy 

ATTN:  POM 

David  Taylor  Naval  Ship  R&D  Ctr. 

ATTN:  Code  1740.5 

ATTN:  Code  L42-3 

ATTN:  Code  2740 

ATTN:  Code  1700,  W.  Murray 

ATTN:  Code  1740,  R.  Short 

ATTN:  Code  177,  E.  Palmer 


31 


DEPARTMENT  OF  THE  NAVY  (Continued) 


DEPARTMENT  OF  THE  NAVY  (Continued) 


Marine  Corp  Dev.  &  Education  Command 

Department  of  the  Navy 

ATTN:  D091,  J.  Hartneady 

Naval  Air  Systems  Command 
ATTN:  F.  Marquardt 

Naval  Construction  Battalion  Center 
ATTN:  Code  L51,  W.  Shaw 
ATTN:  Code  L51,  R.  Odell o 
ATTN:  Code  L51,  S.  Takahashi 
ATTN:  Code  L51,  J.  Crawford 

Naval  Electronic  Systems  Command 
ATTN:  PME  117-21 

Naval  Electronics  Systems  Command 
ATTN:  Commander 

Naval  Explosive  Ord.  Disposal  Fac. 
ATTN:  Code  504,  J.  Petrousky 

Naval  Facilities  Engineering  Command 


Naval  Ocean  Systems  Center 


Office  of  the  Chief  of  Naval  Operations 
ATTN:  OP  981 
ATTN:  OP  982 
ATTN:  OP  982E,  M.  Lenzini 

Strategic  Systems  Project  Office 

Department  of  the  Navy 
ATTN:  NSP-43 

DEPARTMENT  OF  THE  AIR  FORCE 

Aerospace  Defense  Command 

Department  of  the  Air  Force 
ATTN:  XPX 

Air  Force  Armament  Laboratory 
ATTN:  DLYV,  J.  Collins 

Air  Force  Institute  of  Technology 
ATTN:  Library 

Air  Force  Weapons  Laboratory 

Air  Force  Systems  Command 


ATTN:  Code  09M22C 

ATTN: 

NTED 

ATTN: 

SUL 

Material  Command 

ATTN: 

NTES-G,  S.  Melzer 

ATTN:  MAT  08T-22 

ATTN: 

NTE,  M.  Plamondon 

ATTN: 

NTES-C,  R.  Henny 

ATTN 

Code  013,  E.  Cooper 

Ballistic  Missile  Office 

ATTN 

Code  4471 

Air  Force  Systems  Command 

ATTN :  DEB 

Naval 

Postgraduate  School 

ATTN 

Code  0142,  Library 

Ballistic  Missile  Office 

ATTN 

Code  1424,  Library 

Air  Force  Systems  Command 

ATTN:  MMH 

Naval 

Research  Laboratory 

ATTN :  MNNH 

ATTN 

Code  2627 

ATTN 

Code  8440,  G.  O'Hara 

Deputy  Chief  of  Staff 

ATTN 

Code  8403,  R.  Bel  sham 

Research,  Development,  &  Acq 

ATTN 

Code  8440,  F.  Rosenthal 

Department  of  the  Air  Force 

ATTN 

Code  8404,  H.  Pusey 

ATTN:  R.  Steere 

ATTN:  AFRDQSM 

Naval 

Sea  Systems  Command 

ATTN 

SEA-09G53 

Headquarters  Space  Division 

ATTN 

SEA-06 J,  R.  Lane 

Air  Force  Systems  Command 

ATTN 

SEA-033 

ATTN :  DYS 

Naval 

Surface  Weapons  Center 

Headquarters  Space  Division 

ATTN 

Code  U401,  M.  Kleinerman 

Air  Force  Systems  Command 

ATTN 

Code  RIO 

ATTN:  RSS ,  D.  Dowler 

ATTN 

Code  F31 

ATTN 

Code  R14 

Rome  Air  Development  Center 
Air  Force  Systems  Command 

Naval 

Surface  Weapons  Center 

ATTN:  Commander 

ATTN 

Tech.  Library  &  Info.  Services  Branch 

ATTN:  TSLD 

ATTN 

W.  Wi shard 

ATTN:  RBES,  R.  Mair 

Naval 

Weapons  Center 

Strategic  Air  Command 

ATTN 

Code  233 

Department  of  the  Air  Force 

ATTN 

Code  266,  C.  Austin 

ATTN:  XPFS 

Naval  Weapons  Evaluation  Facility 
ATTN :  Code  10 


ATTN:  NR1-STINF0  Library 
DEPARTMENT  OF  ENERGY  CONTRACTORS 


Office  of  Naval  Research 

ATTN:  Code  474,  N.  Perrone 
ATTN:  Code  715 


Los  Alamos  Scientific  Laboratory 
ATTN:  M/S632,  T.  Dowler 
ATTN :  RMS  364 


32 


DEPARTMENT  OF  ENERGY  CONTRACTORS  (Continued) 

Lawrence  Livermore  National  Laboratory 
ATTN:  J.  Thomsen 
ATTN:  M.  Fernandez 
ATTN:  L-200,  T.  Butkovich 
ATTN:  L-90,  R.  Dong 
ATTN:  J.  Goudreau 
ATTN:  L-437,  R.  Schock 
ATTN:  Technical  Information  Dept.  Library 
ATTN:  L-200,  J.  Cortez 
ATTN:  T.  Gold 
ATTN:  L-96,  L.  Woodruff 
ATTN:  L-7,  J.  Kahn 
ATTN:  L-205,  J.  Hearst 
ATTN:  L-90,  D.  Norris 

Oak  Ridge  National  Laboratory 

ATTN:  Central  Research  Library 

Sandia  National  Laboratories 
Livermore  Laboratory 

ATTN:  Library  &  Security  Classification  Div. 

Sandia  National  Laboratories 
ATTN:  A.  Chaban 
ATTN:  W.  Herrmann 
ATTN:  3141 

OTHER  GOVERNMENT  AGENCIES 

Department  of  the  Interior 
Bureau  of  Mines 

ATTN:  Tech.  Lib. 

Federal  Emergency  Management  Agency 

ATTN:  Hazard  Eval.  &  Vul.  Red.  Div. 

U.S.  Nuclear  Regulatory  Commission 

ATTN:  Div.  of  Security  for  L.  Shao 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Acurex  Corp. 

ATTN:  J.  Stockton 

Aerospace  Corp. 

ATTN:  L.  Selzer 

2  cy  ATTN:  Technical  Information  Services 

Agbabian  Associates 

ATTN:  M.  Agbabian 

ATTN:  C.  Bagge 

Analytic  Services,  Inc. 

ATTN:  G.  Hesselbacher 

Applied  Theory,  Inc. 

2  cy  ATTN:  J.  Trulio 

Artec  Associates,  Inc. 

ATTN:  S.  Gill 

AVCO  Research  &  Systems  Group 
ATTN:  Library  A830 
ATTN:  W.  Broding 
ATTN:  J.  Atanasoff 
ATTN:  D.  Henderson 

BDM  Corp. 

ATTN:  R.  Hensley 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
BDM  Corp. 

ATTN:  Corporate  Library 
ATTN:  T.  Neighbors 
ATTN:  A.  Lavagnino 

Bell  Telephone  Labs. 

ATTN:  J.  White 

Boeing  Co. 

ATTN:  Aerospace  Library 

ATTN:  R.  Holmes 

ATTN:  R.  Hager 

ATTN:  M/S  42/37,  R.  Carlson 

ATTN:  J.  Wooster 

ATTN:  R.  Dyrdahl 

Boeing  Co. 

ATTN:  M/S  42/37,  K.  Friddell 

California  Research  &  Technology,  Inc. 

ATTN:  K.  Kreyenhagen 
ATTN:  Library 

California  Research  &  Technology,  Inc. 

ATTN:  D.  Orphal 

Center  for  Planning  &  Rsch.,  Inc. 

ATTN:  R.  Shnider 

Civil  Systems,  Inc. 

ATTN:  J.  Bratton 

University  of  Denver 

ATTN:  Sec.  Officer  for  J.  Wisotski 

EG&G  Washington  Analytical  Services  Center,  Ir 
ATTN:  Library 
ATTN:  Director 

Electric  Power  Research  Institute 
ATTN:  G.  SI  iter 

Electromechanical  Sys.  of  New  Mexico,  Inc. 
ATTN:  R.  Shunk 

Eric  H.  Wang 

Civil  Engineering  Rsch.  Fac. 

ATTN:  D.  Calhoun 
ATTN:  N.  Baum 

Franklin  Institute 

ATTN:  Z.  Zudans 

Gard,  Inc. 

ATTN:  G.  Neidhardt 

General  Dynamics  Corp. 

ATTN:  K.  Anderson 

General  Electric  Co. 

ATTN:  M.  Bortner 

General  Electric  Co. 

ATTN:  A.  Ross 

General  Electric  Company— TEMPO 
ATTN:  DAS  I AC 


33 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

General  Research  Corp. 

ATTN:  B.  Alexander 

H-Tech  Labs,  Inc. 

ATTN:  B.  Hartenbaum 

Honeywel 1 ,  Inc. 

ATTN:  T.  Helvig 

I  IT  Research  Institute 

ATTN:  Documents  Library 
ATTN:  A.  Long i now 

J.  H.  Wiggins  Co.,  Inc. 

ATTN:  J.  Collins 

Kaman  AviDyne 

ATTN:  E.  Criscione 
ATTN:  Library 

Kaman  Sciences  Corp. 

ATTN:  D.  Sachs 
ATTN:  F.  Shelton 
ATTN:  Library 

Karagozian  and  Case 

ATTN:  J.  Karagozian 

Management  Science  Associates 
ATTN:  K.  Kaplan 

Martin  Marietta  Corp. 

ATTN:  G.  Fotieo 
ATTN:  A.  Cowan 

Martin  Marietta  Corp. 

ATTN:  J.  Donathan 

University  of  Massachusetts 
ATTN:  W.  Nash 

McDonnell  Douglas  Corp. 

ATTN:  R.  !lal prin 

Merritt  CASES,  Inc. 

ATTN:  Library 
ATTN:  J.  Merritt 

Mitre  Corp. 

ATTN:  Director 

Nathan  M.  Newmark  Consult.  Eng.  Svcs. 

ATTN:  N.  Newmark 
ATTN:  J.  Haltiwanger 
ATTN:  W.  Hall 

University  of  New  Mexico 

ATTN:  G.  Triandafal idis 

University  of  Oklahoma 
ATTN:  J.  Thompson 

Pacific-Sierra  Research  Corp. 

ATTN:  H.  Brode 

Pacifica  Technology 
ATTN:  G.  Kent 
ATTN:  R.  Bjork 
ATTN:  R.  Allen 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Physics  International  Co. 

ATTN:  F.  Sauer 
ATTN:  R.  Swift 
ATTN:  C.  Vincent 
ATTN:  E.  Moore 
ATTN:  L.  Behrmann 
ATTN:  Technical  Library 

University  of  Pittsburgh 

ATTN:  M.  Wil liras,  Jr. 


R&D  Associates 


P.  Rausch 

Technical  Information  Center 
A.  Field 
R.  Port 
J .  Lewi s 
P.  Haas 


Rand  Corp. 

ATTN:  Library 

Science  Appl ications,  Inc. 

ATTN:  Technical  Library 

Science  Applications,  Inc. 

ATTN:  S.  Os ton 

Science  Applications,  Inc. 

ATTN:  D.  Bernstein 
ATTN:  D.  Maxwell 
ATTN:  R.  Hoffmann 

Science  Applications,  Inc. 

ATTN:  B.  Chambers  III 
ATTN:  W.  Layson 
ATTN:  G.  Binninger 

Southwest  Research  Institute 
ATTN:  A.  Wenzel 
ATTN:  W.  Baker 

SRI  International 

ATTN:  G.  Abrahamson 
ATTN:  W.  Wilkinson 

Systems,  Science  &  Software,  Inc. 
ATTN:  D.  Grine 
ATTN:  Library 

Teledyne  Brown  Engineering 
ATTN:  J.  Ravenscraft 

Terra  Tek,  Inc. 

ATTN:  Library 

Tetra  Tech,  Inc. 

ATTN:  Library 
ATTN:  L.  Hwang 

Texas  ASM  University  System 
ATTN:  H.  Coyle 

TRW  Defense  S  Space  Sys.  Group 
ATTN:  E.  Wong 
ATTN:  F.  Pieper 
ATTN:  G.  Hulcher 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

TRW  Defense  &  Space  Sys.  Group 
ATTN:  A.  Feldman 
ATTN:  P.  Bhutta 
ATTN:  D.  Jortner 

ATTN:  Technical  Information  Center 
ATTN:  B.  Sussholtz 
ATTN:  A.  Narevsky 
ATTN:  P.  Dai 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Weidlinger  Assoc.,  Consulting  Engineers 
ATTN:  J.  McCormick 
ATTN:  M.  Baron 

Westinghouse  Electric  Corp. 

ATTN:  W.  Volz 


Weidl inger  Assoc. ,  Consulting  Engineers 
ATTN:  J.  Isenberg 


35 


